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This paper is devoted to the study of leading twist light cone wave functions of J/^ meson. The 
moments of these wave functions have been calculated within three approaches: potential models, 
nonrelativistic QCD and QCD sum rules. Using the results obtained within these approaches the 
models for the light cone wave functions of leading twist have been proposed. Similarly to the wave 
function of rjc meson the leading twist light cone wave functions of J/ty meson have very interesting 
properties at scales fi > m c : improvement of the accuracy of the model, appearance of relativistic 
tail and violation of nonrelativistic QCD velocity scaling rules. The last two properties are the 
properties of true leading twist light cone wave functions of J/ty meson. 
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I. INTRODUCTION 

Commonly exclusive charmonium production at high energies is studied within nonrelativistic QCD (NRQCD) [1]. 
In the framework of this approach charmonium is considered as a bound state of quark-antiquark pair moving with 
small relative velocity v <C 1. Due to the presence of small parameter v the amplitude of charmonium production can 
be built as an expansion in relative velocity v. 

Thus in the framework NRQCD the amplitude of any process is a series in relative velocity v. Usually, in the most 
of applications of NRQCD, the consideration is restricted by the leading order approximation in relative velocity. 
However, this approximation has two problems which make it unreliable. The first problem is connected with rather 
large value of relative velocity for charmonium: v 2 ~ 0.3, v ~ 0.5. For this value of v 2 one can expect large contribution 
from relativistic corrections in any process. So in any process resummation of relativistic corrections should be done or 
one should prove that resummation of all terms is not crucial. The second problem is connected with QCD radiative 
corrections. The point is that due to the presence of large energy scale Q there appears large logarithmic terms 
(a s log Q/m c ) n , Q 3> m c which can be even more important than relativistic corrections at sufficiently large energy 
( Q ~ 10 GeV). So these terms should also be resummed. In principle, it is possible to resum large logarithms in the 
NRQCD factorization framework [3,3, however such resummation is done rarely. 

The illustration of all mentioned facts is the process of double charmonium production in e + e~ annihilation at 
B-factories, where leading order NRQCD predictions [J, [f| @ are approximately by an order of magnitude less than 
experimental results 0,[J. The calculation of QCD radiative corrections diminished this disagreement but did not 
remove it. Probably the agreement with the experiments can be achieved if, in addition to QCD radiative corrections, 
relativistic corrections will be resummed [To| . 

In addition to NRQCD, hard exclusive processes can be studied in the framework of light cone expansion formalism 
[lllll2| where both problems mentioned above can be solved. Within light cone expansion formalism the amplitude is 
built as an expansion over inverse powers of characteristic energy of the process. Usually this approach is successfully 
applied to excusive production of li ght mesons However recently the application of light cone expansion formalism 
to double charmonium production [13l [Til. [TUlla) allowed one to achieve good agreement with the experiments. 

In the framework of light cone formalism the amplitude of some meson production in any hard process can be 
written as a convolution of the hard part of the process, which can be calculated using perturbative QCD, and process 
independent light cone wave function (LCWF) of this meson that parameterizes nonperturbative effects. From this 
one can conclude that charmonium LCWFs are key ingredients of any hard exclusive process with charmonium 
production. In paper [T^ | the leading twist light cone wave function of ij c meson was studied. This paper is devoted 
to the study of leading twist LCWFs of J/*f? meson. 

The paper is organized as follows. In the next section all definitions needed in our calculation will be given. In 
Section III the moments of LCWFs will be calculated in the framework of Buchmuller-Tye and Cornell potential 
models. Section IV is devoted to the calculation of the moments within NRQCD. QCD sum rules will be applied to 
the calculation of the moments in Section V. Using the results obtained in Sections III-V the models of LCWFs will 
be built in Section VI. In the last section the results of this paper will be summarized. 
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II. DEFINITIONS. 



There are two leading twist light cone wave functions (LCWF) of J /if? meson. The first one is twist 2 LCWF of 
longitudinally polarized J/\f r meson <pt(£i /j). The second one is twist 2 LCWF of transversely polarized J/\P meson 
4>t (£, m) ■ These LCWFs can be defined as follows [l2[ 



(0\Q(z)<T a plz,-z}Q(~z)\J/*(e x=±1 , P ))^ = f T (^)(e aPp - e pPa ) [ d^e^MC,^), (1) 



where the following designations are used: xi , X2 are the parts of momentum of the whole meson carried by quark 
and antiquark correspondingly £ = x\ — x%, p is a momentum of J/^ meson, is an energy scale. The factor [z, —z], 
that makes matrix element ([1]) gauge invariant, is defined as 

[z,-z]=Pexp[ig f dx^A^x)]. (2) 

J — z 

The LCWFs </>l,t(£,m) are normalized as 

J dZ<f> L , T (Z,fx) = l. (3) 

With this normalization condition the constants fo,t are defined as 

(0|g(0) 7Q g(0)|J/*(e A=0 , J3 )} = / L e A=0 , 
(0|Q(0)a^Q(0)|J/*(e A=±1 ,p)) Ai = f T (fJ,)(e a p - e p a ). (4) 

It should be noted here that the local current Q(0)7 Q <5(0) is renormalization group invariant. The local current 
Q(O)cr a0 Q(O) is not invariant. For this reason the constant /t(m) depends on scale fj, as 



^ )= <^) 3b0/ ^' (5) 



but the constant ft does not depend on scale. 

LCWFs <Pl,t{xt^) can be expanded [l2j in Gegenbauer polynomials Cn 2 {£) as follows 



i+ £ ^CmJcS^CO 



(6) 



n=2,4.. 

At leading logarithmic accuracy the coefficients a^ ,T are renormalized multiplicatively 



^• T (m)=(^) £ " /b °^ T (Mo), (7) 



where 

ra+1 



= r 1 



4 £, 



3 V (n+l)(n + 2) 
£ n = J( 4 Ej)' & o = H-^ fl . (8) 

It should be noted here that conformal expansions ([6]) are solution of Bethe-Salpeter equation with one gluon exchange 
kernel [Til ]. 

From equations ©-© it is not difficult to see that at infinitely large energy scale /i — > oo LCWFs 0t,l(£,m) tends 
to the asymptotic form 4> a s(0 = 3/4(1 — £ 2 ). But at energy scales accessible at current experiments the LCWFs 
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4>t.l(€, ar e far from their asymptotic forms. The main goal of this paper is to calculate the LCWFs 0l.t(£, m) °f 
J '/ty meson. These LCWFs will be parameterized by their moments at some scale: 

(£2,t)m = f^CMZ^)- (9) 

It is worth noting that since J/*f? meson has negative charge parity the LCWFs 0l,t(£, A 4 ) are £-even. Thus all odd 
moments y 1 ) equal zero and one needs to calculate only even moments. 
Below the following formulas will be used in our calculation 

(0\Qlu(i^D a ) n Q\J/*(e x=Q ,p)) = f L pA*P) n (®)> 
(0\Q<j^(tz-D a ) n Q\J/^(e x=±1 ,p)) = Mwu-euP M )(zp) n {@), (10) 

where 

D = D-D, D = d -igB a (X a /2). (11) 
These formulas can be obtained if one expands both sides of equations |T|). 

III. THE MOMENTS IN THE FRAMEWORK OF POTENTIAL MODELS. 

In potential models charmonium mesons arc considered as a quark-antiquark system bounded by some static 
potential. These models allow one to understand many properties of chamonium mesons. For instance, the spectrum 
of charmonium family can be well reproduced in the framework of potential models [181 ] . Due to this success one can 
hope that potential models can be applied to the calculation of charmonium equal time wave functions. 

Having equal time wave function of J meson in momentum space ^(k) one can apply Brodsky-Huang-Lepage 
(BHL) [19] procedure and get the LCWFs of leading twist 4>l,t{Z,,p) using the following rule: 

0l,t(C,m) ~ / d 2 fc±Vc(x,k ± ), (12) 

where ?p c (x,\i±) can be obtained from "0(k.) after the substitution 

k ± ^k ± , kz ^( Xl -x 2 )^, M 2 = Mc2+k * . (13) 

Here M c is a quark mass in potential model. In this paper equal time wave function t/>fk) will be calculated in 
the framework of the potential models with Buchmuller-Tye [2(| and Cornell potentials [21(. The parameters of 
Buchmuller-Tye potential model were taken from paper [2(J. For Cornell potential V(r) = —k/r + r/a 2 the calculation 
was carried our with the following set of parameters: k = 0.61, a = 2.38 GeV -1 , M c = 1.84 GeV 22]. 

In paper 17] the moments of leading twist LCWF of r\ c meson were calculated within potential models with these 
potentials. At leading order approximation in relative velocity there is no difference between equal time wave functions 



of r/ c and J/ty mesons. In what follows the moments obtained in paper [17J for the leading twist LCWF of rj c meson 
equal to the moments of LCWFs <Pl,t((, , A 4 ) of J/ty meson. Within this approximation there is no difference between 
<Mf,Ai) and <M£>M)- 

It is worth noting that in paper [2 31 ] the relations between the light cone wave functions and equal time wave 
functions of charmonium mesons in the rest frame were derived. The procedure proposed in paper [23| is similar 
to BHL with the difference: in formula (fT2"]) one must make the substitution d 2 k± — * d 2 k±^/'k 2 + m 2 / '(4m c xiX2)- 
But this substitution was derived at leading order approximation in relative velocity of quark-antiquark motion 
inside the charmonium. At this approximation k 2 ~ 0(v 2 ), Axix 2 ~ 1 + 0(v 2 ) and the substitution amounts to 
d 2 k± — » d 2 k±(l + Q(v 2 )). Thus at leading order approximation applied in [231 ] these two approaches coincide. 

The results of paper [17| are presented in Table I (see this paper for details). In second and third columns the 
moments calculated in the framework of Buchmuller-Tye and Cornell models are presented. It should be noted that 
the moments from Table I were calculated at scale [i ~ 1.5 GeV. It is seen that there is good agreement between these 
two models. 

It should be noted here that the larger the power of the moment the larger contribution form the end point regions 
[x ~ and x ~ 1) it gets. From formulas (|13p one sees that the motion of quark-antiquark pair in this region is 
relativistic and cannot be considered reliably in the framework of potential models. Thus it is not possible to calculate 
higher moments within the potential models. Due to this fact the calculation of the moments has been restricted by 
few first moments. 
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(C) 


Buchmuller-Tye 
model [20.1 


Cornell 
model [21] 


NRQCD 
[24] 


QCD sum rules 

<M£,aO 


QCD sum rules 


n = 2 


0.086 


0.084 


0.075 ±0.011 


0.070 ± 0.007 


0.072 ± 0.007 


n = 4 


0.020 


0.019 


0.010 ±0.003 


0.012 ± 0.002 


0.012 ± 0.002 


n = 6 


0.0066 


0.0066 


0.0017 ±0.0007 


0.0031 ± 0.0008 


0.0033 ± 0.0007 



TABLE I: The moments of LCWFs 4>l{£, a0> 4>t{(,, A 1 ) obtained within different approaches. In the second and third columns 
the moments calculated in the framework of Buchmuller-Tye and Cornell potential models are presented. In the fourth column 
NRQCD predictions for the moments are presented. In last two columns the results obtained within QCD sum rules are shown. 
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IV. THE MOMENTS IN THE FRAMEWORK OF NRQCD. 

In paper [l7| the relations that allow one to connect the moments of leading twist LCWF of rj c meson with NRQCD 
matrix elements were derived 

<e 4 > 

<£ 6 > 

where tp and x + are Pauli spinor fields that annihilate a quark and an antiquark respectively, M c = Mj/^,/2, The 
moments are defined at scale \x ~ M c . 

These relations were derived at leading order approximation in relative velocity. However, as it was noted above 
at this approximation there is no difference between r\ c and J/\& mesons. Moreover, there is no difference between 
LCWFs 4>l (£, /i) and 4>t(£,h)- So the values for the moments of LCWFs ^>l(£, a*), </>t(£> A 4 ) can be taken from paper 

m- 

The results of the calculation of the moments are presented in the fourth column of Table I. The central values of 
the moments and the errors due to the model uncertainty have been calculated according to the approach proposed 
in paper [24| . In addition to the error shown in Table I there is an uncertainty due to higher order v corrections. For 
the second moment one can expect that this error is about 30%. For higher moments this error is larger. 

It is seen from Table I that NRQCD predictions for the second and the fourth moments are in good agreement 
with potential model and there is disagreement for the moment (£ 6 ) between these two approaches. The cause of 
this disagreement is the fact noted above: due to the large contribution of relativistic motion of quark-antiquark pair 
inside quarkonium it is not possible to apply both approaches for higher moments. So one can expect that both 
approaches can be used for the estimation of the values of the second and the fourth moments. The predictions for 
the sixth and higher moments become unreliable. 



V. THE MOMENTS IN THE FRAMEWORK OF QCD SUM RULES. 

In this section QCD sum rules [25l . |2o| will be applied to the calculation of the moments of LCWFs 4>l(£, A*) an d 
0t (£,/•*) [HIH3]- First let us consider LCWF (£,//). To calculate the moments of this LCWF one should consider 
two-point correlator: 

U L (z,q) = i J d 4 xe^(0|TJ (x)J„(0)|0) = (zq) n+2 U L (q 2 ), (15) 

J {x) = Q(x)zQ{x), J„(0) = Q(0)z(iz p D p ) n Q(0), z 2 = 0. 
It is not difficult to obtain sum rules for this correlator (for details see paper fl7j|). 

t^^-^^-PM^+n'zLm, tie) 
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where perturbative and nonperturbative contributions to sum rules ImIIp Cr t(s), H^l Tt (Q 2 ) can be written as 

n^ rt (Q 2 ) = n( m) (g 2 ) + n( m) (g 2 ) + n( m) (Q 2 ), (is) 



24^ 7-1 V 4 7 (4m 2 + Q 2 (l-£ 2 )) m+2 

n< m) (Q 2 ) = -%^!l m 2 ( m 2 + 3m+2) r^ra+s^ (i-e 2 r +1 



(4m 2 + Q 2 (l-C 2 ))' 



n ( m '(Q 2 ) = %^>(„a_„)( m +i) /' ^T~ 2 (i-e 2 )"^ 



384^ ( 4TO 2 + Q2 (1 _ C 2)y m + 2 - 

Here Q 2 = — q 2 , m c and (a s G 2 ) are parameters of QCD sum rules. 

To calculate the moments of LCWF 0t(£,aO one should consider the correlator: 

U T (z,q) = i f d*xe^{0\TJp{x)J£(p)\0) = (zg)" +2 n T (g 2 ), (19) 



J^x) = Q(x)(a^z»)Q(x), J£(0) = Q(0){a^z u )(i Z PD p ) n Q(0), z 2 = 0. 
The sum rules for this correlator can be written as 



frM = I r rf ImII pert (s) (m) 

(m 2 J/3 , + Q2 )m+ 1 „ {s + Q2 )m+ 1 + npcrtV 



^ t^zu +k;uq% (so) 



where perturbative and nonperturbative contributions to sum rules (|20p are given by expressions (|17[) . (|18p except 
that the expression for 11^ (Q 2 ) should be replaced by 



n w (Q2) = <^> (m+ i) (_ r+ ^) r - 2( i-a 



(1 - £2xm+2 

1 g J 3,. (21) 



(4m 2 + Q 2 (l-£ 2 )) 



In the original paper [26| the method QCD sum rules was applied at Q 2 — 0. However, as it was shown in paper 
28], there is large contribution of higher dimensional operators at Q 2 = which grows rapidly with m. To suppress 
this contribution sum rules (fT6| . (|20| will be applied at Q 2 = Am 2 . 

In the numerical analysis of QCD sum rules the values of parameters m c and (a s G 2 /ir) will be taken from paper 

M- 

m c = 1.24 ± 0.02 GeV, (— G 2 ) = 0.012 ± 30% GeV 4 . (22) 

First sum rules ([TBI) . (12"0)) will be applied to the calculation of the constants f 2 T . It is not difficult to express the 
constants f 2 T from equations (|16p . (f2"0"|) at n = as functions of m. For too small values of m (m < mi) there are 
large contributions from higher resonances and continuum which spoil sum rules (|16p. (|20p . Although for m 3> mi 
these contributions are strongly suppressed, it is not possible to apply sum rules for too large m (m > m 2 ) since the 
contribution arising from higher dimensional vacuum condensates rapidly grows with m and invalidates approximation 
of this paper. If mi < m 2 there is some region of applicability of sum rules (|16|) . (|20|) [mi, m 2 ] where both resonance 
and higher dimensional vacuum condensates contributions are not too large. Within this region f 2 T as a functions 
of m vary slowly and one can determine the values of these constants. Applying approach described above one gets 

fl = 0.170 ± 0.002 ± 0.004 ± 0.016 GeV 2 , 

/| = 0.167 ± 0.002 ± 0.003 ± 0.016 GeV 2 . (23) 

The first error in (|23|) corresponds the variation of the constants / 2 T within the region of stability. The second and 
the third errors in (|23p correspond to the variation of the gluon condensate (a s G 2 ) and the mass m c within ranges 
From the results (|2"3"|) one sees that the main errors in determination of the constants f 2 T result from the 
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variation of the parameter m c . This fact represents well known property: high sensitivity of QCD sum rules to the 
mass parameter m c . 

Next let us consider the second moments T ) in the framework of QCD sum rules. One way to find the values of 
(£i t) i s ^0 determine the values of j\ T ) from sum rules (fT6|) . ([20]) at n = 2 and then extract (£f T ). However, 
as it was noted above, this approach suffers from high sensitivity of right side of equations (fT6|) , (|20|) to the variation 
of the parameter m c . Moreover, the quantities f 2 T ) include not only the errors in determination of T ), but 
also the errors in f 2 T ,. To remove these disadvantages the ratios of sum rules at n = 2 and n = 0: f 2 T (£f, t)I ft t 
will be considered. The moments T ), (££ T ) will be considered analogously. Applying standard procedure one gets 
the moments of LCWF </>l(£, ^): 

= 0.070 ± 0.002 ± 0.007 ±0.002, (24) 
= 0.012 ±0.001 ± 0.002 ±0.001, 
(d) = 0.0031 ±0.0002 ±0.0008 ±0.0002, 

and moments of LCWF 0t(£, A 4 ) 

(£t) = 0.072 ±0.002 ±0.007 ± 0.002, (25) 
(£t) = 0.012 ±0.001 ± 0.002 ±0.001, 
(Ct> = 0.0033 ±0.0002 ±0.0007 ±0.0003. 

The first error in (|24|) . (f2~5|) corresponds the variation within the region of stability. The second and the third errors 
in (f24|) . (|25|) correspond to the variation of the gluon condensate (a s G 2 ) and the mass m c within ranges ([22]) . It is 
seen that, as one expected, the sensitivity of the ratios f 2 t)//E t t° the variation of m c is rather low. The main 
source of uncertainty is the variation of gluon condensate (a s G 2 ). In the fifth and sixth columns of Table I results 
([2"5|) are presented. The errors in Table I correspond to the main source of uncertainty — the variation of gluon 
condensate (a s G 2 ). 

In the calculations of the correlators (fT5|) . fH]) characteristic virtuality of quark is ~ (Am 2 + Q 2 )/m ~ m 2 . So the 
values of the moments (|2l)|) are defined at scale ~ m 2 . 

From Table I it is seen that the larger the number of the moment n the larger the uncertainty due to the variation 
of vacuum gluon condensate. This property is a consequence of the fact that the role of power corrections in the sum 
rules (|16p . (|20p grows with n. From this one can conclude that there is considerable nonperturbative contribution to 
the moments (££ T ) with large n what means that nonperturbative effects are very important in relativistic motion 
of quark-antiquark pair inside the meson. The second important contribution to QCD sum rules (|16p . (|20|1 at large n 
is QCD radiative corrections to perturbative part n por t((5 2 ). Unfortunately today one does not know the expression 
for these corrections and for this reason they are not included to sum sules ([175)) . ([2Tj|) . One can only say that these 
corrections grow with n and, probably, the size of radiative corrections to the ratios /f T (££ T ) / 1 f 2 T is not too big 
for not too large n. Thus one can expect that QCD radiative corrections will not change dramatically the results for 
the moments n = 2 and n — 4. But the radiative corrections to (££ T ) may be important. 

It is interesting to compare the moments of leading twist LCWF of r\ c meson calculated in paper [TtJ 

(£2j = 0.070 ± 0.002 ±0.007 ±0.003, (26) 
(^j = 0.012 ±0.001 ± 0.002 ±0.001, 

_ 0.0032 ±0.0002 ±0.0009 ±0.0003, 

with the results of this section. It is seen that there is no significant difference between the moments of leading 
twist LCWF of r] c meson and the moments of LCWFs 4>l{£,, /■*)> <Ar(C) A*)- The difference is within the error of the 
calculation. 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 



Fig.l: The LCWF <j> L (29]) at scales fi = 1.2 GeV,/n = 10 GeV,/t 2 = 100 GeV,/i 3 = oo. 

VI. THE MODEL FOR THE LCWFS OF J/* MESON. 

In paper pjj the model of leading twist LCWF of rj c meson was proposed: 

#£,M= Mo) = c(/3)(l -ftexp^-JL-y (27) 

where c(/3) is a normalization coefficient, the constant j3 = 3.8 ± 0.7, the scale /io = 1.2 GeV. This function allows one 
to reproduce the results (|26[) with rather good accuracy. The moments of this wave function are 

(£ 2 ) = 0.070 ±0.007, (28) 
(£ 4 ) = 0.012 ±0.002, 
(£ 6 ) = 0.0030 ±0.0009. 

At central value j3 = 3.8 the constant c((3) ~ 62. 

As it was noted in the previous section the accuracy of the calculation does not allow one to distinguish leading 
twist LCWF of rj c meson from LCWFs of J r /\f r meson. For this reason the model ([27jl will be used for the LCWFs 
4>l(L M = Mo), <Mf , M = Mo) 

<Mf,M= Mo) = <M£,M = Mo) = </>(£, M = Mo)- (29) 

In this expression the functions /z = /io), 0t(£,M = Mo) are defined at scale fi = /io- It is not difficult 

to calculate these functions at any scale /i > /io using conformal expansions ([6]). The LCWFs </>l(£,m) a t scales 
/io = 1.2 GeV, /ii = 10 GeV, /i2 = 100 GeV, /i3 = oo are shown in Fig. 1. The moments of LCWFs </>l(£, /i) at scales 
/io = 1.2 GeV, /ii = 10 GeV, n% — 100 GeV, /i3 — oo are presented in second, third, fourth and fifth columns of Table 
II. The plot and the moments of LCWF <^(£, /i) will not be shown since this function practically does not deviate 
from 

As was noted in paper [17| model ((29|) has some interesting properties. For instance, let us consider the LCWF 
</>l(£, /i) (similar consideration can be done for </V(£, m))- From conformal expansion ([6]) one can derive the expressions 
that determine the evolution of the moments: 





i 

5 4 


hi v!2 

« 2 (m)^, 




(30) 


= 


3 

35 
1 

21 


+ 4(m)|- 

£/ \ 1 2 


r , , 120 r. , 64 

4VW 1001 b "^2145 





Similar relations can be found for any moment. Further let us consider the expression for the second moment 
In this paper the value has been found with some error at scale /i = /io- This means that the value of the 
coefficient a L (m = Mo) was found with some error. The coefficient decreases as scale increases. So the error in 
a\ and consequently in decreases as scale increases. At infinitely large scale there is no error in at all. 
The calculations show that the error 10% in at scale /i = /io decreases to 4% at scale /i = 10 GeV. Applying 
relations (|30p it is not difficult to show that similar improvement of the accuracy happens for higher moments. The 
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improvement of the accuracy allows one to expect that model (j29|) at larger scales will be rather good even if QCD 
radiative corrections to results (j2"S"|) are large. 

From Fig. 1 one sees that LCWF at scale /i = /io practically vanishes in the regions 0.75 < |£| < 1. In this region 
the motion of quark-antiquark pair is relativistic and vanishing of LCWF in this region means that at scale /i = /io 
charmonium can be considered as a nonrelativistic bound state of quark-antiquark pair with characteristic velocity 
v 2 ~ 1/(3 ~ 0.3. Further let us regard the function 4>l(£,,h = /io) as a conformal expansion ([6]). To get considerable 
suppression of the LCWF in the region 0.75 < |£| < 1 one should require fine tuning of the coefficients of conformal 
expansion a^(fi = /io)- The evolution of the constants (especially with large n) near fi = /io is rather rapid (see 
formulas ([7]) and (|8|)) and if there is fine tuning of the constants at scale /i = /io this fine tuning will be rapidly broken 
at larger scales. This property is well seen in Table II and Fig. 1. From Fig. 1 it is seen that there is relativistic tail 
in the region 0.75 < |£| < 1 for scales \i = 10, 100 GcV which is absent at scale /i = /i . Evidently this tail cannot 
be regarded in the framework of NRQCD. This means that, strictly speaking, at some scale charmonium can not be 
considered as nonrelativistic particle and the application of NRQCD to the production of charmonium at large scales 
may lead to large error. Although in the above arguments the model of LCWF (|29|) was used it is not difficult to 
understand that the main conclusion is model independent. 

According to the velocity scaling rule [l[ the moments of LCWF depend on relative velocity as ~ v n . It is not 
difficult to show that the moments of LCWF (|29|) satisfy these rules. Now let us consider the expressions that allows 
one to connect the coefficients of conformal expansion with the moments (££). These expressions for the moments 
(££)> (££) are given by formulas (|30|) . It causes no difficulties to find similar expressions for any moment. From 
expressions (f3"0j) one sees that to get velocity scaling rules: (££) ~ v n at some scale one should require fine tuning of 
the coefficients at this scale. But, as was already noted above, if there is fine tuning of the coefficients at some 
scale this fine tuning will be broken at larger scales. From this one can conclude that velocity scaling rules are broken 
at large scales. 

Consider the moments of LCWF ([29]) at infinite scale. It is not difficult to find that 

3 

(Cl,t)m=oo = t — — Y7 — —zz- (31) 
(n + l)(n + 3) 

From last equation one can find that (££. T ) does not scale as v n as velocity scale rules [lj require. Thus scaling rules 
obtained in paper [l[ are broken for asymptotic function. Actually one does not need to set the scale /t to infinity to 
break these rules. For any scale /i > /io there is a number uq for which the moments (£™ T ), n > uq violate velocity 
scaling rules. This property is a consequence of the following fact: beginning from some n = n the contribution of 
the relativistic tail of LCWF, that appears at scales /t > /io, to the moments becomes considerable. 

The amplitude T of any hard process with charmonium meson production can be written as a convolution of LCWF 
$(£) with hard kernel H(£) of the process. If one expands this kernel over £ and substitute this expansion to the 
amplitude T one gets the results: 

T = f dCH(Om = E ^T^(D- (32) 

If one takes the scale /i ~ /io in formula (|52")l , than moments (£") scale according to the velocity scaling rules ~ v n 
and one gets usual NRQCD expansion of the amplitude. However due to the presence of the scale of the hard 
process \ih ^ Mo there appear large logarithms log/t/j//io which spoil NRQCD expansion (|3"2"|) . To remove this large 
logarithms one should take fi ~ /i/j. But at large scales velocity scaling rules are broken and application of NRQCD 
is questionable. 

In papers [HI, [23|, [2i| it was proposed different models of LCWFs of J/\& and rj c mesons. It is interesting to compare 
the models proposed in these papers with model (|2"!)j) . Such comparison was done in paper [ijj and it will not be 
discussed in this paper. 

VII. CONCLUSION. 

In this paper the moments of leading twist light cone wave functions (LCWF) of J/ meson have been calculated 
within three approaches. In the first approach Buchmullcr-Tye and Cornell potential models were applied to the 
calculation of the moments of LCWFs. In the second approach the moments of LCWFs were calculated in the 
framework of NRQCD. In the third approach the method QCD sum rules was applied to the calculation of the 
moments. The results obtained within these three approaches are in good agreement with each other for the second 
moment (£ 2 ). There is a little disagreement between the predictions for the fourth moment (£ 4 ). The disagreement 
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<o 


<H£, /xo = 1.2 GeV) 


<?!>(£, Mi = 10 GeV) 


4>{£,H2 = 100 GeV) 


<K£, A*3 = oo) 


n = 2 


0.070 


0.12 


0.14 


0.20 


n = 4 


0.012 


0.040 


0.052 


0.086 


n = 6 


0.0032 


0.019 


0.026 


0.048 



TABLE II: The moments of LCWF proposed in this paper at scales fi = 1.2 GeV,A*i = 10 GeV,fj, 2 = 100 GeV,^ 3 = oo 
are presented in second, third, fourth and fifth columns. 



between the approaches becomes dramatic for the sixth moment (£ 6 ). The cause of this disagreement consists in the 
considerable contribution of relativistic motion of quark- antiquark pair inside J/^Sf meson to higher moments which 
cannot be regarded reliably in the framework of potential models and NRQCD. The approach based on QCD sum 
rules is more reliable, especially for higher moments since it does not consider J/'F-meson as a nonrelativistic object. 
The main problem of QCD sum rules is that since there is no expressions of radiative corrections to sum rules one 
does not know the size these corrections. However, one can expect that QCD radiative corrections will not change the 
results for the moments n = 2 and n = 4 dramatically. As to the sixth moment, the contribution the QCD radiative 
corrections in this case may be important. 

The moments of leading twist LCWFs of J/\I/ meson have been compared with the moments of LCWF of r/ c . It 
was found no significant difference between the moments of leading twist LCWF of r\ c meson and the moments of 
LCWFs 4>l(£,, a0> aO- The difference is within the error of the calculation. For this reason the model of LCWF 
of rj c meson was taken as a model for leading twist LCWFs </>l(£, a0> </>t(£, m) °f J/^ meson. As it was shown in paper 
[l7| this model has some interesting properties: 

1. Due to the evolution © the accuracy of the moments obtained within model (f2"9")l improves as the scale rises. 
For instance, if the error in determination of the moment T ) is 10% at scale // = A*o = 1-2 GeV, at scale a* = 10 
GeV the error is 4%. For higher moments the improvement of the accuracy is even better and there is no error at all 
at infinite scale (i = oo. The improvement of the accuracy allows one to expect that model (f2"9"]l will be rather good 
even after inclusion of the QCD radiative corrections. 

2. At scale [i ~ /Xo the LCWFs can be considered as wave functions of nonrelativistic object with characteristic 
width ~ v ~ 0.3. Due to the evolution, at larger scales relativistic tail appears. This tail cannot be considered in 
the framework of NRQCD and, strictly speaking, at these scales J/^> meson is not a nonrelativistic object. 

3. It was found that due to the presence of high momentum tail in the LCWFs at scales fi > /io there is violation 
of velocity scaling rules obtained in paper [lj . More exactly, for any scale \i > [io there is a number no for which the 
moments (££ T ),n > uq violate NRQCD velocity scaling rules. 

Actually, the last two properties are properties of real LCWFs of J/^ meson. 

The author thanks A.V. Luchinsky and A.K. Likhoded for useful discussion and help in preparing this paper. This 
work is partially supported by Russian Foundation of Basic Research under grant 04-02-17530, Russian Education 
Ministry grant RNP-2.2.2. 3.6646, CRDF grant MO-011-0, president grant MK-1863.2005.02 and Dynasty foundation. 
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